Reconstruction of condensed magnetoexciton droplet in a trap in strong magnetic 

fields 
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We investigate theoretically the Bose-Einstein condensation of trapped magnetoexcitons in a two- 
layered system with one layer containing electrons and the other layer containing holes. We have 
studied the spatial variations of the condensate density in the droplet of electrons and holes. We 
find that the shape of the electron and hole densities may change due to the competition between 
repulsive electron-electron/hole-hole interaction and confinement potential of the trap. Our mean 
field calculations show that when the confinement strength is strong enough the condensate density 
is peaked at the edge of the droplet, and as the confinement strength weakens the condensate density 
displays one inner peak and another peak at the edge. For much weaker confinement potential the 
condensate density may display one broad peak. 

PACS numbers: 78.55.Cr, 03.75.nt, 73.20.Mf, 78.67.De 



I. INTRODUCTION 

The Bose-Einstein condensation (BEC) of excitons is 
more likely to occur when (a) electrons and holes are 
spatially separated, (b) a strong magnetic field is present, 
and (c) bosons are confined in a trap. When electrons and 
holes are spatially separated exciton lifetime increases. 
The lifetime depends on the wavcfunction overlap be- 
tween electrons and holes. In coupled quantum wells the 
wavefunction overlap can be controlled by an external 
electric field applied perpendicular to the quantum wells. 
It was shown theoretically that in such a system exci- 
tonic condensate can be realized in strong magnetic fields 
nearly at all filling factors when the separation between 
the electron and hole layers is less than the magnetic 
length 1.) . Only for large separations fractional quantum 
Hall states or Wigner crystal states may be realized. Dif- 
ferent aspects of excitonic BEC in two dimensional sys- 
tems in strong magnetic fields have been investigated over 
the years [H 11,11,11,011,011,11, liolllll • As the observa- 
tion of atomic BEC demonstrates [lj the conditions for 
BEC are improved if bosons are confined in a trap. Sev- 
eral groups have recently investigated excitonic BEC in 
a trap in the absence of a magnetic field. Experimentally 
the possibility of excitonic BEC was explored in traps of 
double quantum wells 0, Q, 0] . Theoretically exciton 
BEC in a trap was investigated in a mean field theory in- 
cluding electron-electron and hole-hole interactions p^. 
Signatures of BEC in angular distribution of photolumi- 
nescence have been also explored theoretically 17]. 

In this paper we investigate theoretically the BEC of 
excitons when both a strong magnetic field and a trap are 
present in double quantum wells. The magnetic field is 
applied perpendicular to the layers and the confinement 
potential traps particles laterally, i. e. in the layers. We 
consider the strong magnetic field limit where only the 
states in the lowest Landau level are relevant. The trap 



potential will split the degenerate lowest Landau level 
states. What is interesting about strong magnetic field 
limit is that the particle density cannot exceed a certain 
value l/2ir£ 2 , where £ is the magnetic length. In bulk 2D 
this density is achieved when the lowest Landau level if 
completely filled, i. e. when the Landau level filling factor 
is one. A strong confinement potential pushes particles 
in each layer to the center of the potential and the par- 
ticle density takes the largest possible value except near 
the edge of such a droplet. Similar effect is not present in 
the absence of a strong magnetic field. We will call this 
state the maximum density droplet (MDD). When the 
confinement potential becomes weaker the particle den- 
sity can take a smaller value at some place in the droplet. 
Such a reconstruction for MDD 0, E, [2(j in electron 
single dots was investigated both theoretically and ex- 
perimentally by several groups in strong magnetic fields 

m e e mi m m m m m 113. The purpose f 

this paper is to investigate how a reconstruction of the 
MDD affects the condensed magnetoexcitons. 

We have explored this problem within a BCS theory. 
In the MDD condensed magnetoexcitons are found near 
the edge of the droplet, see Fig. [21 We find that as 
the strength of the confinement potential weakens this 
MDD becomes unstable and density depletion starts to 
occur at the interior of the droplet. In this state con- 
densed magnetoexcitons are found both in the interior of 
the droplet and near the edge of the droplet, see Fig|3] 
When the confinement potential is weakened further the 
density depletion of electrons and holes increases. This 
behavior is similar to how the electron MDD reconstructs 
in electron single dots j2jj. However, we find that, unlike 
the electron single dot case, where the number of elec- 
trons depleted increases discontinuously from one, two, 
and etc., [2p| . the density depletion increases continu- 
ously with the decrease in the strength of the confine- 
ment potential. When density changes from the uniform 
state are more severe the order parameter may display 
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one broad peak. 

This paper is organized as follows. In Sec. II we de- 
scribe our model and its Hamiltonian. A gap equation 
for the condensate is derived in Sec. III. In Sec. IV we show 
solutions of this gap equation. Conclusions are given in 
Sec.V. 



II. MODEL 



We will derive a gap equation for BEC of two- 
dimensional magnetoexcitons in a trap described 
above. Here we are particularly interested in the com- 
petition between the confinement potential and the re- 
pulsive particle-particle interactions. We thus include 
electron-electron and hole-hole interactions in the Hamil- 
tonian explicitly. Our system then consists of pseudo 
spin-up and -down particles in two dimensions under a 
strong magnetic field. The Hamiltonian of the system is: 



An artificial trap for excitons in double quantum wells 
can be set up by changing locally the well width on one 
side of double quantum well [T(j- An in-plane random 
potential may also provide a local potential minimum 
which can host electrons and holes [TH ITEj. A parabolic 
potential may be created by applying inhomogeneous 
stress Another possibility of a trap is self-assembled 
quantum dots 28] . 

We take a simple model for the confinement potential 
of such traps and consider electrons and holes to be con- 
fined in their respective layers by parabolic quantum dot 
potentials, described by a parameter fi: ^m e fl 2 r 2 and 
^m^f^r 2 , where m e and are the electron and hole 
masses. A magnetic field B is applied perpendicular to 
the 2D layers. (Hereafter the perpendicular direction will 
be taken to be the z-axis) . We confine our attention here 
to the strong magnetic field limit, where f2/w C)(T < 1, 
with cJ C)ff = eB /m a c and a = e,h. We assume that the 
magnetic field is so strong that both electron and hole 
systems are spin-polarized. In this limit the symmetric 
gauge single-particle eigenstates are conveniently classi- 
fied by a Landau level index n and an angular momentum 
index m. Here we consider the strong field limit so n = 
and m = 0, 1, 2, ... Electrons with angular momentum m 
have a wavefunction 



<Pm(r) 



1 
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(1) 



while holes with angular momentum —m have a wave- 
function 



<f>- m (r) = 



20mm! \i 



(2) 



Here z = x + iy and £ — ^/hc/eB. Note that these wave- 
functions are independent of electron and hole masses. 
An electron with angular momentum m has a wave- 



function 



(r) 



rotating counter clockwise. A 



hole with angular momentum — m has a wavefunction 
0_ TO (r) ~ e lrn< ^ rotating clockwise. To stress the similar- 
ity between our approach and the treatment of supercon- 
ductivity we will rename electrons and holes as pseudo 
spin-up and -down particles. The single particle orbitals 
in the n = level have energies t a m = huj c _ a /2+'y a (m+l), 
where j a = m a £l 2 i 2 , uj Cj(T — eB/m a c, and a =|, | is the 
pseudospin index. The term 7 ct (to + 1) represents the 
splitting of the degenerate states of the lowest Landau 
level by the potential of the trap. The electron and hole 
creation operators are c^* and c_ m >. 



(3) 



where 



H ? = H e '» c ™T Cm T ( 4 ) 

II I 

rrii ,m' 

= ^el m cL mi c_ mi (5) 

m 

+ \ X! Vsi-m'^-m^-mi 



-m 2 i 



are respectively the Hamiltonian for the electrons and 
the holes, and 

H c -h = -- ^ V d (m' 1 ,-m 2l m 1 ,-m 2 ) 

XC li' 1 ] C -rn' 2 i C - m 2l c mih ( 6 ) 

Hh-c = -- Y Vd{-m' 1 ,m' 2 ,-m l7 m 2 ) 

xc -m' 1 i c m;,T Cm2 T c - m ii (^) 
describe the electron- hole interactions. Here J2 m m' 
stands for X)m x m 2 m' m' an d Vd is defined to be positive 
so we put minus signs in front of 1/2 in 7i e -h and Ti.h- e - 
The matrix elements for intra dot electron-electron/hole- 
holc and inter-dot electron-hole interactions are, respec- 
tively, 

V s (m' 1 ,m' 2 ,mi,m 2 ) = /^ r i J d 2 r 2 V s {r! ~ r 2 ) 

x C;( r i)Ci( r 2)0mi(ri)<Wr2) (8) 

and 

Vd(m' 1 ,m' 2 ,mi,m2) = J d 2 r r J d 2 r 2 V d (ri - r 2 ) 

xCi (ra)^mi (ri)0m 2 (r a ). (9) 

Note that V s (r) is just the Coulomb interaction and 



V d (f) 



= , where d is the inter layer dis- 



e^Jr 2j rd 2 

tance and e is the dielectric constant of the semicon- 
ductor. One can show that V^m'i, — m' 2 , mi, — m 2 ) = 
Vd(-m' 2 , m'x, -m 2 , mi). 
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III. GAP EQUATIONS AND SELF ENERGIES 



where 



In a bulk 2D when the filling factor v ^ 1 the excitonic 
condensate is described well by the mean field BCS state 



-hf\ 



(m)+e l HF (m) - 2fi+. 



I*> = f[(um + v m cl n ^_ ml )\0). 



(10) 



This state represents a uniform density condensate with 
constant the particle occupation number, — v. Small 
system exact diagonalization studies indicate the over- 
lap between this mean field uniform BCS and exact 
groundstate wave functions is close to one for small val- 
ues of d/i. Also a trial wavefunction leads to identical 
physical results^ as mean field field results^. 

We apply the same trial wavefunction to describe the 
excitonic BEC in a trap 163 • The expectation values of 
the different terms in the Hamiltonian with respect to 
| \&) may be conveniently evaluated by using 

c - m l = u m a_ m ^ — v m a m ^ 

(11) 

The a's satisfy anticomutation relations. The BCS state 
1$) represents the vacuum state: a m j|<I') = and 
a m|.W — 0- Using the usual Hartree-Fock (HF) de- 
coupling we find the expectation values 

(H e ) = £{4 + -£[F s (m,m',m,TO') 



V s (m, to', to', m)]v^ n ,}v 2 



(12) 



(Kb) = J2i e L + 2 [V s (-m,-m,-m,-m') 

711 — 77l f —0 

- V s (—m, —to', —to', —m)]v^ n ,}v'^ n (13) 

{H e -h} + (Hh-c) = 

- X! V d (m,~m,m / ,~m / )vl l vl l , 

771,771' =0 
771c 

- X! ^(to, -to, m',-m')u m v m u m rv m '. 



m, m'— 



(14) 



We can derive gap equations by minimizing the total en- 
ergy (H c ) + (Hb) + (W e -h) + (Wh-c):- 



A m = - Vd (to, -to, to', -to') 



(15) 




(16) 

(17) 
(18) 



Here we define the chemical potential fi + = (/ij + //j.)/2. 
In solving the gap equation we must determine /i + self- 
consistently. The occupation number of single particle 
states is v^. The HF effective single particle energy of a 
particle with pseudo spin a is 

e° HF [m) = e° m + E^m) + E x (m) + E^m), (19) 

where is the single particle confinement potential en- 
ergy and 

Efl- iS (m) = V s (to, to', to, Tn')v m , , 

Ex(to) = - ^ y s (to, to', to', to)^, , 

^H,d( m ) = -/;Vd(m, -m, to', -m^v^,. (20) 

The physical meaning of each term is as follows. Con- 
sider, for example, the HF effective single particle en- 
ergy of a pseudo spin-up particle. The first and second 
terms then represent the Hartree and exchange self en- 
ergy corrections due to the presence of other pseudo spin- 
up particles. The third term represents the Hartree self 
energy correction due to the presence of pseudo spin- 
down particles. Note that when the inter layer distance 
d = then E x (to) = E^.^to) since V s (m,m' ,to' ,to) = 
Vd{m, —to, to', —to'). 

IV. RESULTS 

We present numerical solutions of the gap equations, 
Eqs. H15fl - I|16|) . f° r a trap with the electron layer and 
hole layer separated in the z-direction by a distance d = 
0.5£. For this inter layer distance BEC is expected in 
bulk 2D systems for filling factors less than one Q. The 
spread of the wave functions of electrons and holes in 
the z-direction is assumed to be negligibly small. We 
investigate the ideal case with T = and equal electron 
and hole masses m e = rrih — 0.067too. In our study we 
have equal number of electrons and holes, N e = Nh = N. 

A. Test of the mean field theory 

Before we present the numerical results we test the 
accuracy of our HF approach. In the limit where the 
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FIG. 1: Results for the occupation numbers of single particle 
states n m in the absence of electron-hole interaction. Cir- 
cles are HF results and diamonds are exact result (Solid and 
dotted lines are guides to the eye). Here hQ — 3meV and 
N = 20. 
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FIG. 2: Left panel: Particle occupation numbers v m (solid 
line) are plotted as a function of quantum number m, which 
label single particles states. u^ n — I — are also plotted 
(dotted line). Right panel: The order parameter g(r) is plot- 
ted as a function of r. Here Ml — 2.35meV and N = 30. For 
" 2 and g(r) for d — are also shown (dashed 



comparison v„ 
lines) . 



electron-hole interaction is absent electrons and holes 
are independent of each other and exact diagonaliza- 
tion results are known |2jJ . In Fig. ^ we have com- 
pared HF occupation numbers of single particle states 
with those of exact results (In the HF approach occupa- 
tion numbers are just n m = vf n ). The following param- 
eters are used: N — 20, Ml — 3meV, fi + = 5meV, and 
B = 5.9T. We see that when the effect of quantum fluc- 
tuations are included the suppression in the occupation 
numbers become more broader and less deeper as a func- 
tion of to. Also the positions of the minimum of the den- 
sity depletion are somewhat different between the two. 
Nonetheless we see that HF correctly captures qualita- 
tively the right physics of reconstruction near the MDD. 
The total amount of electron or hole density depletion in 
the interior of the droplet is exactly one in both cases: 
V 15 u 2 - 1 



B. Order parameter of reconstructed droplet 

The order parameter, i. e., the condensate density, is 

g{r) = (* T (r)* x (r)) 

= ^2 ( c m-[C- m >i)(j) m (r)(j)- m '(r), (21) 
mm' 

where * T (r) = E c mT^™( r ) and *l( r ) = 
Y)c— m \<p— m (r) are electron and hole field opera- 
tors, res pec tively. Since the optical selection rule of 
excitons |29j requires to = to' we have 



ffOO = ^2(c„ l TC^ ml )\(j) m (r)\ 1 



(22) 



Since 



we have 



^ Crn^C— m[ — ^m^m — o ' (^^) 

4 EJm 



](r) = 22v m u m \(j> m (r)\' 1 



(24) 



Note that |</>m(r)| 2 is peaked at r ~ \/2(m + 1)£ with a 
width £ (see EqEJ. 

Before we show the results for the order parameter let 
us first investigate when the MDD is realized. Fig. [5] 
displays the calculated and u 2 m for Ml = 2.35TOel^, 
N = 30, d = 0.U, and B = 5T. The chemical potential 
is /i + = b.&2meV. The electron and hole occupation 
numbers, v^, are one except near the edge region to « 
to c = 29. Note that the radius of the edge of the droplet 
is r c = y / 2(m c + 1)1. The electron or hole density of this 
MDD is given by 

n e , h (r) =Y,v 2 m \&n(r)\ 2 . (25) 

This density looks approximately like a step function 



n e ,h(r) 




(26) 



Thus, in a parabolic trap with equal number of elec- 
trons and holes a uniform density state can be realized in 
the strong magnetic field limit with the particle density 
1/2tt<? 2 . 

Let us investigate the condensate order parameter for 
this state. Near the edge, where to c ~ 29, both and 




FIG. 3: Same as in Fig. [5] but with a smaller value of the FIG. 5: Same as in Fig. |5|but with Ml = 2.0meV. 

strength of the confinement potential hfl = 2.27meV . 




u 2 m are non-zero, and, consequently, the order parameter 
v m u m is non-zero, see Fig. A rough shape of g{r) may 
be deduced as follows. The product v m u m is roughly 
given by a delta function 6 m>nic . It then follows from 
Eq. i|24|) that g(r) ~ \4 > mA r )\ 2 > which is peaked near the 
edge r c . In this MDD electrons and holes are surrounded 
by condensed excitons near the edge of the droplet Q . We 
have also calculated and g(r) for an increased value for 
the strength of electron-hole interaction, corresponding 
to d = (see dashed lines). We notice that the width of 
the ring where the order parameter is non-zero becomes 
larger compared to the case of d = 0.5£, see the dashed 
curve in the right panel of Fig. [3 

We now show the results on how the uniform MDD 
reconstructs as the strength of the confinement potential 
decreases. We will call these reconstructed states nearly 
uniform droplets. When a smaller value Ml = 2.21meV 
is used vL near m — deviates from one and wL starts to 



take non zero values, see Fig. [3]. The self-consistent value 
of the chemical potential is /i + = 5.76meV. In this case 
the order parameter develops two peaks. This feature can 
be understood by investigating the product of v m u m : it is 
non zero near m = and m — 29, and two peaks are ex- 
pected. The quantity J^m u m = 0-92 gives the number of 
electrons or holes added from the interior of the droplet 
to the edge. As Kl is reduced further to 2.2meV the 
value of is almost zero at m = and J2 m u m = 1-85 
(see Fig. 0]). This means that more particles are added 
from the interior of the droplet to the edge. The value 
of the chemical potential is /i + = 5.63meV. When Ml is 
reduced even further to 2.0meV the occupation number 
is almost zero, not at m = 0, but, at finite value of 
m (see Fig. 0). In this case the value of the chemical 
potential is ^ + = 5.35meV. Note that near m = 3 the 
occupation numbers is nearly zero while v? m takes a 
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FIG. 7: The effective mean field single particle energies as 
functions of to. Solid line is for hQ. = 2.35meV ( see Fig. 
and the dashed line is for hQ — 2.0meV ( see Fig. |^J. For 
both cases N = 30. Horizontal lines represent the values of 
the chemical potentials. 



maximum value. The sum is u m = 4.23. We con- 
clude from results in Figs. [21 - G3 that with decreasing 
value of Ml the total density added from the interior of 
the droplet to the edge increases continuously. Similar 
trend also exists in the case of electron single dots, but 
there it changes discontinuously and takes integer values. 
A non-uniform state is also investigated for Ml — 2meV, 
fi+ = 3.85meV, and N = 20. Fig. displays the cal- 
culated v^, and the order parameter for this case. 
The occupation numbers are non-uniform with the 
average value 1/2 (For larger values of Ml studied above 
the average occupation number was one) . The order pa- 
rameter is maximum between the center and the edge of 
the droplet, but does not exhibit a double peak as be- 
fore. We believe that the results of the HF approach are 
less reliable in describing non uniform states than nearly 
uniform states. The results of Fig[f)]should thus be taken 
as rather approximate results. 



C. Mean field single particle energies 

A strong confinement potential pushes particles in each 
layer to the center of the potential and the particle den- 
sity of the droplet takes the largest possible value except 
near the edge of such a MDD. When the confinement po- 
tential becomes weaker some particles move from inside 
of the droplet to the edge, i.e., a reconstruction takes 
place, and, consequently, the size of the droplet expands 
in the plane of 2D system. It is not simple to predict 
exactly where this density depletion takes place in the 
droplet since it is a consequence of a non-trivial inter- 
play between the confinement potential, Hartree poten- 
tial, exchange potential, and electron-hole interactions. 



m 



FIG. 8: Various terms of the mean field single particle energy 
for HQ = 2.27meV and N = 30 are plotted (The correspond- 
ing order parameter is plotted in Fig. . Top panel: Solid 
fine represents the intra layer Hartree self energy, Ejj |S (to), 
while dashed line represents the confinement single particle 
energy. Middle panel: Solid line represents the intra layer ex- 
change self energy, Ex (to), while dashed line represents the 
inter layer Hartree self energy Y*H,d( m )- Bottom panel: Solid 
line represents the effective mean field single particle energy, 
e a HF {m), as functions of to while dashed line represents the 
chemical potential. 
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FIG. 9: Same as in Fig. |H1 but for the non-uniform state. 
The corresponding order parameter is shown in Fig. |S| 



The mean field single particle energy eJ^ F (m), given 
in Eq. I|19l) . reflects this interplay to some degree. A 
rough shape of may be obtained from the shape of 
the mean field single particle energy e a HF {m): In the ab- 
sence of electron-hole pairing correlation single particle 
states with e^ F (m) smaller than the chemical potential 
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are occupied with probability one. However, since we 
do have pairing in the BEC this is only approximately 
correct. The mean field single particle energies e^ F (m) 
are plotted for Ml — 2.35meV and 2.0meV in Fig. [7| 
The corresponding are shown in Fig [3 and Fig [SJ re- 
spectively. As the strength of the confinement potential 
decreases the local maximum of e a HF (to) moves from the 
center of the droplet and the width of this local maximum 
becomes broader. The resulting occupation numbers 
are approximately consistent with the shape of e a HF (m). 
Even for a non-uniform droplet the position of the peak 
in t'jjp (to) coincides with the position of the suppression 
of v 2 ^. see the peak near to = 7 in e t j IF {m) (FigJSJ and 
the dip in near m — 7 (Fig. EJ). 

The various self energies appearing in e^ F (m) are plot- 
ted for Ml — 2.27meV in Fig. [S] The presence of a peak 
in e a HF {rn) near to = is related to the peaks in the 
self energies Ex(to) and E^ ^to). The corresponding 
occupation number decreases near m = , see Fig. 
H3 (We notice that T,H,d{m) is almost flat in the region 
6 < to < 20. £# j£ i(TO) displays qualitatively the same 
dependence as the exchange self energy Hx(m)- When 
d = they are actually identical and are flat for to < m c 

E3). 

V. CONCLUSIONS AND DISCUSSION 

We have investigated the shape of the condensed mag- 
netoexcitons in circular traps in strong magnetic fields. 
Our model consists of a parabolic confinement in lateral 
directions and a delta function like confinement along the 
perpendicular axis. We have applied a mean field theory, 
which is expected to be qualitatively correct. We find 
that as the strength of the confinement potential weak- 
ens, or equivalently as B increases, the uniform state 
becomes unstable and density depletion starts to occur 
in the interior of the droplet. We found that the amount 
of density depletion increases continuously with the de- 
crease in the strength of the confinement potential. As a 
consequence of these reconstructions the order parameter 
changes from displaying one peak at the edge to display- 
ing one inner peak and another peak at the edge for de- 
creasing confinement strength. When density depletions 
are more severe, ie, when the confinement potential is 



rather weak, the order parameter may display one broad 
peak. 

The reconstruction may be observed experimentally 
since the spatial shape of the order parameter changes. 
The structure in the order parameter g(r) may be inves- 
tigated experimentally by observing angular distribution 
of photolumincscence since it reflects the Fourier trans- 
formed order parameter in the momentum space. An ex- 
pression for how the PL angular profile depends on the 
order parameter is given explicitly in Keeling et al |l7j . 
They point out that it may provide a diagnostic for the 
existence of BEC. 

Note that decreasing Ml and increasing B has the same 
effect. Experimentally either of these two parameters 
may be varied. This can be understood as follows. The 
strength of the confinement strength enters as a dimen- 
sionless parameter in units of the Coulomb energy scale: 
7 CT = r y a /(e 2 /e£). The reconstruction of the shape of the 
droplet will depend on this parameter, which is a func- 
tion of B and Ml: jcr ~ {Ml) 2 / B z / 2 . It may be easier 
experimentally to change B than Ml. 

Also it is desirable to increase the particle numbers one 
order of magnitude and investigate the shape of the con- 
densate droplet. However, this is a non-trivial task since 
the matrix elements of the particle-particle interactions 
are difficult to calculate numerically when the quantum 
number to are large. 

We comment on effects that are not investigated in 
this paper. It would be interesting to study how the cor- 
rections to the Hartree Fock theory affect quantitative 
aspects of the shape of the order parameter. This can 
be investigated by performing numerical exact diagonal- 
ization. Finite quantum well width corrections of the 
electron and hole wavefunctions will weaken intralayer 
interactions and hence favor excitonic states over FQHE 
states. Realistic valence band structure effects [2^ can 
be included through the matrix elements of the intra and 
inter layer Coulomb interaction. Corrections to the low- 
est Landau level approximation may have quantitative 
relevance [30|. 

This work is supported by Korea Research Founda- 
tion grant KRF-2003-015-C00223, and by grant No.(R01- 
1999-00018) from the interdisciplinary research program 
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